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ABSTRACT 

Turbine  blade  ehdwall  heat  transfer  measurements 
are  given  for  a  range  of  Reynolds  and  Mach  numbers. 
Data  were  obtained  for  Reynolds  numbers  based  on  in¬ 
let  conditions  of  0.5  and  1.0  x  10^,  for  isentropic  exit 
Mach  numbers  of  1.0  and  1.3,  and  for  freestream  tur¬ 
bulence  intensities  of  0.25%  and  7.0%.  Tests  were  con¬ 
ducted  in  a  linear  cascade  at  the  NASA  Lewis  Transonic 
Turbine  Blade  Cascade  Facility.  The  test  article  was  a 
turbine  rotor  with  136®  of  turning  and  an  axial  chord 
of  12.7  cm.  The  large  scale  allowed  for  very  detailed 
measurements  of  both  flow  field  and  surface  phenom¬ 
ena.  The  intent  of  the  work  is  to  provide  benchmark 
quality  data  for  CFD  code  and  model  verification.  The 
flow  field  in  the  cascade  is  highly  three-dimensional  as 
a  result  of  thick  boundary  layers  at  the  test  section  in¬ 
let.  Endwall  heat  transfer  data  were  obtained  using  a 
steady-state  liquid  crystal  technique. 

LIST  OF  SYMBOLS 

A  -  heater  area  [cm^] 

Cp  -  specific  heat  [J /kg*K] 

Cx  -  blade  axial  chord  [cm] 

Ec  -  Eckert  number,  Ec  =  U^/{CpAT) 
h  -  heat  transfer  coefficient  [W/m^-K] 
k  -  thermal  conductivity  [W/m-K] 
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electrical  current  [Amps] 
heater  length  in  y-direction  [cm] 
longitudinal  integral  turbulence  length 
scale  [cm] 

cascade  mass  flow  rate  [kg/s] 

Mach  number 
blade  pitch  [cm] 
pressure  [Pa] 

Prandtl  number 
heat  transfer  rate  [Watts] 
recovery  factor,  r  =  Pr^/^ 
heater  resistance  [Ohms] 
gas  constant  [J/kg-K] 

Reynolds  number,  Rec^  =  pUinCx/p 
blade  span  [cm] 

Stanton  number 
thickness  [cm] 
temperature  [K] 

Turbulence  intensity 
total  velocity  [m/s] 
electrical  voltage  [Volts] 
chordwise  (axial)  direction 
pitchwise  (tangential)  direction 
spanwise  (radial)  direction 
pitch  angle  (x-y)  plane  [deg.] 
specific  heat  ratio,  7  =  1.4 
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99%  boundary  layer  thickness  [cm] 
emissivity 

dynamic  viscosity  [kg/s-m] 
density  [kg/m^ 

Stefan-Boltzmann  constant 

adiabatic  wall  temperature 

bulk  heat  addition 

liquid  crystal 

exit  freestream  value 

inlet  freestream  value 

isentropic  value 

heater 

evaluated  at  Eckert’s  reference 

temperature 

surface 

total  conditions 


INTRODUCTION 

As  the  dependence  on  computational  analyses  for 
giero dynamics  and  heat  transfer  increases  in  the  design 
of  turbomachinery,  so  does  the  need  for  detailed,  bench¬ 
mark  quality  experimental  data.  These  data  are  re¬ 
quired  for  validation  of  existing  turbomachinery  CFD 
analyses,  and  to  provide  data  for  improved  approaches 
to  modeling  the  flow  field  characteristics.  Because  of 
the  sensitivity  of  the  life  of  high  pressure  turbine  ro¬ 
tors  to  metal  temperature,  validation  of  the  procedures 
used  to  predict  these  temperatures  is  essential.  Vali¬ 
dation  of  CFD  analyses  should  be  done  at  conditions 
which  approach  engine  operating  conditions  as  closely 
as  possible.  Several  experiments  have  obtained  de¬ 
tailed  data  for  turbine  rotor  geometries.  Langston  et 
al.  (1977),  Gregory-Smith  and  Graves  (1983),  Mar- 
chal  and  Sieverding  (1977),  Yamamoto  (1987a,b),  and 
Gregory-Smith  et  al.  (1988)  presented  detailed  results 
of  aerodynamic  measurements  for  rotor  blades,  tested 
in  linear  cascades  at  relatively  low  Mach  numbers,  Mee 
et  al.  (1992a)  and  Mee  et  al.  (1992b)  presented  mea¬ 
surements  for  a  rotor  linear  cascade  tested  at  transonic 
conditions  obtained  in  a  short  duration  blowdown  fa¬ 
cility.  Heat  transfer  measurements  for  turbine  rotor 
geometries  were  obtained  by  Graziani  et  al.  (1980), 
Goldstein  and  Spores  (1988)  in  linear  cascades,  and  by 
Blair  (1994)  in  a  rotating  cascade.  All  of  the  heat  trans¬ 
fer  experiments  were  done  in  large  scale  test  facilities, 
but  were  at  low  Mach  number.  To  increase  confidence 
in  the  ability  of  computational  analyses  to  predict  tur¬ 
bine  rotor  heat  transfer  under  actual  engine  conditions, 
there  exits  a  need  to  obtain  both  aerodynamic  and  heat 
transfer  data  at  high  Mach  number  conditions. 


A  Transonic  Turbine  Blade  Cascade  Facility  was 
built  at  the  NASA-Lewis  Research  Center  to  provide 
detailed  aerodynamic  and  heat  transfer  measurements 
at  high  Mach  numbers.  The  facility  is  a  linear  cascade. 
It  was  recognized  that  a  linear  cascade  cannot  simulate 
rotation  effects,  and  that  the  chosen  design  does  not  al¬ 
low  for  tip  clearance  effects.  Nonetheless,  a  linear  cas¬ 
cade  was  chosen  over  an  annular  cascade  for  a  number 
of  reasons.  The  primary  reason  was  to  obtain  detailed 
measurements  in  a  large  scale  facility  at  transonic  flow 
conditions.  The  physical  size  of  the  blading  was  de¬ 
termined  by  the  amount  of  air  that  could  be  moved  at 
transonic  speed  through  the  cascade.  An  annular  cas¬ 
cade  of  engine  typical  hub-to-shroud  and  aspect  ratios 
would  require  several  times  more  flow  rate  than  a  lin¬ 
ear  cascade  for  blades  of  the  same  chord  length.  Also, 
a  linecLT  cascade  provides  easier  access  for  instrumenta¬ 
tion  than  an  annular  cascade.  From  a  computational 
aspect,  a  linear  geometry  minimizes  the  effects  of  the 
three-dimensional  mesh  on  the  CFD  solution,  thus  giv¬ 
ing  a  clearer  analysis  of  the  actual  flow  solver. 

Some  noteworthy  features  of  the  linear  cascade  fa¬ 
cility  are:  its  ability  to  achieve  transonic  Mach  num¬ 
bers,  its  high  blade  count  to  ensure  good  periodicity, 
and  its  large  scale.  The  large  scale  facilitates  highly 
localized  and  detailed  flow  and  heat  transfer  measure¬ 
ments.  Because  of  the  large  scale  {Cx  =  12.7  cm),  and 
transonic  flow  conditions,  inlet  Reynolds  numbers  up 
to  one  million  can  be  achieved  with  an  inlet  pressure 
near  one  atmosphere.  The  purpose  of  the  research  in 
this  facility  is  to  obtain  a  detailed  aerodynamic  and 
heat  transfer  database  for  Reynolds  and  Mach  numbers 
relevant  to  high  pressure  turbine  applications.  The  re¬ 
sults  of  the  investigation  are  available  for  CFD  code  and 
model  validation. 

The  rotor  tested  in  the  cascade  is  of  constant  cross 
section,  so  that  the  geometry  of  the  cascade  is  two- 
dimensional.  It  was  desired  to  provide  verification  data 
for  three-dimensional  CFD  analyses.  The  three  dimen¬ 
sionality  of  the  flow  field  is  achieved  by  not  bleeding 
off  the  endwall  boundary  layers.  Giel  et  al.  (1996) 
show  that  when  the  endwall  boundary  layers  are  not 
bled  off,  the  spanwise  variation  of  surface  pressure  on 
the  rotor  was  of  the  same  order  as  the  spanwise  varia¬ 
tion  in  surface  pressure  experienced  in  an  actual  engine. 
For  some  tests  a  square  bar  grid  was  placed  upstream 
of  the  blade  row  to  generate  high  levels  of  freestream 
turbulence.  Without  the  turbulence  grid  present,  the 
measured  inlet  freestream  turbulence  was  about  0.25%. 
With  the  grid  installed,  the  freestream  turbulence  was 
in  excess  of  7%. 
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Endwall  heat  transfer  data  was  obtained  using  a 
steady-state  liquid  crystal  technique.  The  liquid  crys¬ 
tal  heat  transfer  measurement  technique  was  chosen  be¬ 
cause  of  the  good  spatial  resolution  that  it  offers.  Sev¬ 
eral  versions  of  the  technique  have  been  used  for  gas  tur¬ 
bine  related  measurements.  A  transient  technique  was 
used  by  Martinez-Botas  et  al.  (1994)  while  a  steady- 
state  technique  was  used  by  Hippensteele  and  Russell 
(1988). 

The  objective  of  the  current  paper  is  to  present  the 
results  of  endwall  heat  transfer  measurements  that  were 
made  in  the  Transonic  Turbine  Blade  Cascade.  Data 
were  obtained  for  eight  different  flow  conditions.  The 
results  from  the  eight  cases  will  be  presented  to  illus¬ 
trate  the  effects  of  varying  Reynolds  number,  exit  Mach 
number,  and  the  inlet  turbulence  grid.  Heat  transfer 
measurements  were  made  at  inlet  Reynolds  numbers  of 
0.5  and  1.0  x  10®.  Tests  were  made  at  pressure  ratios 
corresponding  to  exit  isentropic  Mach  numbers  of  1.0 
and  1.3.  Tests  were  conducted  at  high  and  low  inlet 
turbulence  levels. 


DESCRIPTION  OF  FACILITY 

The  initial  description  of  the  facility  was  given  by 
Verhoff  et  al.  (1992).  Subsequently,  the  cascade  inlet 
section  was  analyzed  and  redesigned  (Giel  et  al.  1994) 
to  improve  the  inlet  flow  pitchwise  uniformity.  Figure  1 
shows  an  overall  view  of  the  facility,  including  the  re¬ 
designed  forward  inlet  boards. 


High  pressure  air  at  ambient  temperature  is  throt¬ 
tled  to  give  a  maximum  total  pressure  of  approximately 
one  atmosphere  at  the  inlet  to  the  test  section.  The  high 
Reynolds  number  tests  had  an  inlet  total  pressure  near 
one  atmosphere,  and  the  low  Reynolds  number  tests 
had  an  inlet  total  pressure  near  one-half  atmosphere. 
The  air  was  discharged  into  an  exhaust  header  main¬ 
tained  at  a  nominal  pressure  of  15.9  kPa  (2.3  psia).  A 
valve  between  the  test  section  and  exhaust  header  was 
used  to  throttle  the  flow  to  establish  the  desired  exit 
Mach  number.  As  seen  in  Fig.  1,  the  test  section  is 
mounted  on  a  large  disk.  This  disk  can  be  rotated  to 
give  a  range  of  incidences.  All  of  the  results  described 
were  obtained  at  the  design  inlet  flow  angle  of  63.6  de¬ 
grees.  This  inlet  angle  gave  a  flow  turning  of  about 
136  degrees.  Upstream  inlet  boards  were  installed  as 
described  in  Giel  et  al.  (1994),  but  to  prevent  shock  re¬ 
flections  in  the  downstream  section,  no  exit  tailboards 
were  used.  A  highly  three-dimensional  flow  field  was 
obtained  in  the  blade  passages  by  allowing  boundary 
layers  to  develop  in  the  long  inlet  section  upstream  of 
the  cascade.  Aerodynamic  probe  data  and  blade  load¬ 
ing  data  both  verify  the  existence  of  strongly  three- 
dimensional  passage  flow. 

The  blade  shape  used  in  these  tests,  while  generic 
in  nature,  was  designed  to  simulate  the  flow  character¬ 
istics  of  a  high  specific  work  rotor.  Pertinent  details  of 
the  blade  and  of  the  cascade  are  given  in  Table  1.  A 
plan  view  of  the  test  section  is  shown  in  Fig.  2.  The 
figure  shows  the  five  0.025  mm  (0.001  inch)  thick  In¬ 
conel  heater  sheets  and  their  positions  relative  to  the 
blades.  Passage  5  is  the  primary  measurement  passage. 
The  heated  section  begins  upstream  of  the  blade  lead¬ 
ing  edge  plane  at  — 0.3Cx  and  extends  downstream  to 
2.1  Cx-  Each  of  the  heaters  are  L  =  2.1  Cx  wide  in  the 
^-direction.  The  figure  also  shows  the  outline  of  the 
6.35  cm  (2.50  in.)  thick  plexiglas  window  that  was  used 
to  gain  optical  access  to  the  test  section.  The  window 
formed  the  endwall  opposite  the  instrumented  endwall. 
The  camera  fields  of  view  are  also  shown  in  the  figure. 
Camera  1  was  used  primarily  to  obtain  data  in  the  blade 
leading  edge  region  while  Camera  2  was  used  primarily 
in  the  passage  and  trailing  edge  regions.  Camera  3  had 
a  wider  angle  lens  which  allowed  a  view  of  the  entire 
test  section.  It  was  thus  used  to  obtain  data  in  the  far 
downstream  region,  as  well  as  in  the  upstream  region 
near  the  heater  leading  edge.  Camera  I’s  view  of  this 
region  was  shadowed  by  the  frame  holding  the  plexigleis 
window.  All  of  the  cameras  were  located  2.0  m  from  the 
instrumented  endwall.  Cameras  1  and  2  had  135  mm 
focal  length  lenses  while  Camera  3  had  an  85  mm  focal 
length  lens. 


Fig.  1  Overall  view  of  Transonic  Turbine 
Blade  Cascade  test  section 
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Table  1  Blade  and  cascade  parameters 
and  dimensions 


Geometric  parameter 

Value 

axial  chord 

pitch 

span 

true  chord 
stagger  angle 
throat  diameter 
throat  area:  1  passage 
leading  edge  diameter 
trailing  edge  diameter 
turbulence  grid 

12.70  cm  (5.000  inches) 
13.00  cm  (5.119  inches) 
15.24  cm  (6.000  inches) 
18.42  cm  (7.250  inches) 
41.54" 

3.358  cm  (1.393  inches) 
53.94  cm^  (  8.360  in^) 
2.657  cm  (1.046  inches) 
0.518  cm  (0.204  inches) 
2.54  cm  square  bar 

Flow  parameter 

Value 

Inlet  Rec^ 

Exit  iZec* 

Inlet  Mis 

Exit  Mis 

Inlet  5/  -  no  grid 
with  grid 
Inlet  flow  angle 

Design  flow  turning 

1.0  X  10®  or  0.5  X  10® 

1.8  X  10®  or  0.9  X  10® 
0.38 

1.32  or  0.97 

3.2  cm  (1.2  inch) 

2.0  cm  (0.8  inch) 

63.6* 

136® 

The  outer  wall  of  the  instrumented  test  section 
was  3.8  cm  (1.5  in.),  thick  stainless  steel  to  withstand 
a  pressure  difference  of  one  atmosphere.  There  was 
2.3  cm  (0.90  in.)  of  low  thermal  conductivity  {kfoam  = 
0.029  W/mK),  high  temperature,  closed  cell  foam  be¬ 
tween  the  outer  wall  and  the  heaters.  The  foam  was 
coated  with  0.3  cm  (0.10  in.)  of  gelcoat  epoxy  to  pro¬ 
vide  a  smooth  surface.  ‘T’-shaped  electrical  bus  bars 
extended  through  the  endwall  assembly.  They  were 
flush  with  the  endwall  and  passed  through  the  outer 
wall.  The  inner  ends  of  the  bus  bars  are  evident  in 
Fig.  2  at  both  ends  of  all  five  heaters.  The  Inconel  foil 
was  tack  welded  to  the  bus  bars  along  these  ends.  The 
bus  bars  were  thermally  as  well  as  electrically  conduc¬ 
tive,  and  were  therefore  placed  outside  the  test  pas¬ 
sage.  The  rest  of  the  foil  was  bonded  to  the  gelcoat 
surface  using  double-sided  adhesive  tape.  The  heaters 
were  connected  in  series  outside  the  test  section. 

The  endwall  heater  assembly  was  first  painted  with 
flat  black  lacquer  paint  and  then  a  uniform  grid  of  white 
dots  were  painted  on  with  a  template.  The  dots  were 
used  in  the  data  reduction  process  to  determine  spatial 
location.  The  locations  of  the  reference  dots  are  shown 
in  Fig.  2. 

Because  the  blades  were  made  of  stainless  steel, 
those  in  contact  with  the  heaters  needed  to  have  electri¬ 
cally  insulated  ends.  RTV  silicone  was  chosen  to  coat 
the  ends  of  the  blades.  A  thickness  of  0.4  mm  (1/64  in.) 


was  determined  from  a  heat  conduction  analysis  to  best 
balance  the  insulating  effect  of  the  blade,  which  would 
tend  to  allow  overheating,  and  the  fin  cooling  effect  of 
the  blade,  which  would  tend  to  overcool  the  endwall 
region  near  the  blade. 

The  optional  74%  open  area  square  bar  turbulence 
grid  is  also  shown  in  Fig.  2.  The  grid  consisted  of  one 
square  25  mm  x  25  mm  (1  in.  x  1  in.)  pitchwise  bar 
that  extended  between  the  inlet  boards  at  midspan  and 
three  square  spanwise  bars  located  on  150  mm  centers. 

MEASUREMENT  TECHNIQUES 
Inlet  Flow  Measurements  As  shown  in  Fig.  2,  aero¬ 
dynamic  probe  measurements  were  made  on  a  plane  one 
axial  chord  upstream  of  the  blade  leading  edge  plane. 
The  measurements  covered  three  blade  pitches  and  ex¬ 
tended  from  near  the  endwall,  zfs  =  8.3  x  lO""*,  to  just 
above  midspan,  z/s  —  0.54.  Details  of  the  measure¬ 
ment  techniques  are  given  by  Giel  et  al.  (1996).  Time 
mean  flow  measurements  were  made  with  a  calibrated 
3-hole  boundary-layer  probe  near  the  endwall  and  a  cal¬ 
ibrated  5-hole  pitch-yaw  probe  away  from  the  endwall. 
Turbulence  measurements  were  made  with  a  constant 
temperature  hot  wire  anemometer. 

Endwall  Static  Pressure  Measurements  The  local 
isentropic  Mach  number  is  needed  to  calculate  the 
endwall  surface  Stanton  number  distribution.  Surface 
static  pressure  data  was  obtained  with  a  separate  end- 
wall,  instrumented  only  with  pressure  taps.  The  test 
passage  and  two  adjacent  passages  were  each  instru¬ 
mented  with  87  pressure  taps  per  passage.  The  taps 
extended  axially  from  x/Cx  =  —0.4  to  x/Cx  =  2.1, 
and  were  arranged  in  21  rows  with  either  4  or  5  taps 
per  row.  Pressure  measurement  taps  were  placed  in  all 
passages,  but  away  from  the  test  passage  fewer  mea¬ 
surements  were  made.  The  locations  of  the  taps  will  be 
shown  along  with  the  data  in  a  following  section.  More 
details  of  the  endwall  static  pressure  measurement  tech¬ 
niques  are  given  by  Giel  et  al.  (1996). 

Heat  Transfer  Measurements  The  steady-state  liq¬ 
uid  crystal  technique  was  chosen  for  the  present  study 
for  several  reasons.  Primarily,  the  method  is  rela¬ 
tively  insensitive  to  substrate  non-homogeneity  (e.g.  air 
bubbles,  thickness  variations,  etc.).  Furthermore,  the 
isotherms  can  be  photographed  with  higher  quality  still 
photographs  rather  than  with  video  that  has  poorer 
spatial  and  color  resolution.  The  following  is  a  brief 
overview  of  the  technique:  Thin  Inconel  sheets  were 
placed  on  the  endwall  so  as  to  cover  the  surface  of  the 
test  blade  passage.  A  suspension  of  liquid  crystals  was 
painted  on  the  Inconel  sheets.  Electrical  power  was 
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supplied  to  the  Inconel  sheets  to  generate  a  spatially 
uniform  heat  flux.  For  each  liquid  crystal,  there  is  a 
very  narrow  temperature  band  at  which  a  well-defined 
color  is  visible.  Where  the  color  is  observed,  the  sur¬ 
face  temperature  is  known.  Since  the  electrical  power 
is  known,  the  heat  transfer  coefficients  can  be  deter¬ 
mined.  Varying  the  electric  power  changes  the  location 
of  the  liquid  crystal  isotherm.  When  steady-state  con¬ 
ditions  were  achieved,  the  locations  of  the  isotherms 
were  photographically  recorded.  Finally,  the  locations 
of  the  isotherms,  along  with  endwall  isentropic  Mach 
number  distributions,  were  used  to  generate  contours 
of  Stanton  number.  Details  of  all  of  the  steps  in  this 
process  will  now  be  discussed. 

Two  micro-encapsulated  chiral  nematic  liquid  crys¬ 
tals  with  yellow-band  temperatures  of  Tjc,i  =  37.9±2®C 
(100.2  ±  3.6"F)  and  Tic,2  =  64.5  ±  3"C  (148.1  ±  5.4"F) 
were  mixed  and  sprayed  onto  the  surface.  The  crys¬ 
tal  temperatures  were  calibrated  using  the  surface- 
mounted  thermocouples  shown  in  Fig.  2.  As  described 
by  Moffat  (1990),  micro-encapsulated  crystals  tend  to 
be  less  sensitive  to  illuminating  and  viewing  angles  than 
non-encapsulated  crystals.  While  the  optical  axes  of 
all  three  cameras  were  orthogonal  to  the  planar  end- 
wall,  strobe  lights  used  for  illumination  were  mounted 
at  various  angles.  To  best  illuminate  the  views,  two 


photographs  were  taken  with  each  of  the  three  cameras, 
one  with  a  strobe  to  the  right  of  the  camera  and  one 
with  a  strobe  to  the  left.  No  differences  in  isotherm  lo¬ 
cation  were  ever  detected  between  the  two  photographs. 

The  voltage  drop,  V,  across  each  of  the  five  heaters 
was  measured  independently.  The  current,  7,  through 
the  entire  circuit  was  measured  with  three  redundant 
0.1%  accuracy  shunts.  From  these  voltage  and  cur¬ 
rent  measurements,  the  power,  qhtr  =  V/,  and  the 
resistance,  Tl  =  V/7,  were  calculated,  displayed,  and 
recorded. 

One  disadvantage  of  the  steady-state  liquid  crys¬ 
tal  method  is  that  longer  run  times  are  required  for 
thermal  equilibration  to  be  achieved.  After  changing 
the  heater  power  level,  the  steady-state  condition  was 
examined  by  monitoring  the  isotherm  position  on  a  re¬ 
mote  video  display.  The  resistance  of  the  heaters  was 
also  monitored.  Inconel  was  chosen  as  the  heater  mate¬ 
rial  because  of  its  low  coefficient  of  thermal  resistivity. 
The  signals  giving  the  heater  voltage  and  current  had 
high  signal-to-noise  ratios.  Even  though  the  calculated 
resistance  varied  by  less  than  0.5%  over  the  entire  range 
of  heater  power  settings,  variations  in  resistance  with 
temperature  could  be  monitored.  Stability  of  the  re¬ 
sistance  was  also  used  as  an  indicator  of  when  steady 
state  conditions  were  achieved. 
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When  steady-state  conditions  were  achieved, 
35  mm  color  slide  photographs  were  taken  of  the  end- 
wall  surface  with  three  cameras.  Each  camera  was  posi¬ 
tioned  to  view  a  subset  of  the  entire  surface,  with  some 
overlap  between  camera  views  as  described  previously 
and  as  shown  in  Fig.  2.  The  end  wall  was  illuminated 
using  high  speed  strobe  lights.  The  short  duration  of 
the  strobe  light  prevented  radiative  heating  of  the  liq¬ 
uid  crystals.  Two  strobes  were  used  for  each  camera, 
and  were  positioned  to  prevent  blind  spots  caused  by 
reflections  or  shadows.  Pressures,  temperatures,  and 
other  data  were  recorded  simultaneously  with  the  pho¬ 
tographic  record.  The  resulting  photographic  slides 
were  then  projected  onto  paper  marked  with  the  dot 
pattern  shown  in  Fig.  2.  The  dot  pattern  in  the  slide 
image  was  aligned  with  the  dot  pattern  on  the  paper. 
Composite  drawings  of  the  yellow  line  positions  were 
made.  Each  camera  and  crystal  temperature  had  a 
separate  drawing.  Flow  condition  data,  time,  crystal 
temperature,  and  power  level  were  associated  with  each 
line.  These  lines  were  then  digitized.  The  coordinates 
and  their  associated  data  were  stored  for  further  pro¬ 
cessing.  Typically,  about  5000  digitized  points  were 
stored  for  each  flow  condition  run.  A  sample  of  the 
digitized  point  distribution  is  shown  in  Fig.  3.  The 
full  range  of  surface  Stanton  numbers  was  typically  ob¬ 
tained  from  25  different  power  levels.  Typical  power 
settings  ranged  from  200  to  1200  Watts  in  increments 
of  15  to  100  Watts.  The  particular  power  settings  used 
depended  on  the  inlet  ambient  temperature. 


Fig.  3  Sample  digitized  data  point  distribution 

The  coordinates  of  each  digitized  point,  along  with 
the  data  associated  with  the  point,  were  read  into  a  data 
reduction  program.  The  purpose  of  the  program  was  to 
calculate  Stanton  numbers  at  each  digitized  point,  but 
this  calculation  required  knowledge  of  the  local  adia¬ 
batic  wall  temperature.  The  first  step  in  the  procedure 
was  thus  to  read  in  CFD  calculations  of  endwall  sur¬ 
face  static  pressure,  which  will  later  be  shown  to  agree 
well  with  the  experimental  data.  The  calculated  val¬ 
ues  provided  significantly  better  spatial  resolution  than 
the  experimental  pressure  distribution.  Isentropic  sur¬ 


face  Mach  numbers  were  determined  at  each  digitized 
point  by  interpolation  from  the  CFD  calculations.  Af¬ 
ter  determining  M/5,  Taw  and  thus  the  Stanton  number 
could  be  determined  at  each  of  the  digitized  points.  The 
Stanton  and  Mach  numbers  were  then  interpolated  onto 
a  uniform,  rectangular  grid  to  generate  contour  plots. 
For  all  eight  flow  conditions,  a  176  x  41  grid  was  used  on 
each  heater.  Three  point  linear  interpolation  was  used, 
and  care  was  taken  never  to  interpolate  across  the  blade 
boundaries.  Since  the  Mach  numbers  originated  as  cal¬ 
culated  values,  they  had  extremely  fine  spatial  resolu¬ 
tion.  Contour  plots  of  the  interpolated  Mach  numbers 
were  also  created  for  verification.  Comparing  contour 
plots  of  the  original,  input  Mach  number  distribution 
with  contour  plots  generated  after  the  Mach  numbers 
had  been  first  interpolated  onto  each  digitized  point, 
and  then  interpolated  onto  the  rectangular  grid,  vali¬ 
dated  the  interpolation  procedure.  If  the  two  contour 
plots  did  not  agree,  it  was  determined  that  this  was  a 
region  of  insufficient  data  in  which  no  Stanton  number 
contours  should  be  shown.  For  all  of  the  test  cases  this 
situation  arose  only  near  the  edges  of  the  measurement 
area. 

The  range  of  electric  power  supplied  to  the  end- 
wall  heaters  was  limited  so  that  they  became  neither 
too  hot  nor  too  cool.  The  heaters  were  assembled  in 
a  highly  stretched  condition  at  room  temperature.  At 
low  temperatures  the  thin  heaters  could  stretch  beyond 
their  elastic  limit,  and  then  wrinkle  when  warmed.  The 
heaters  could  become  too  cool  at  low  power  levels,  es¬ 
pecially  near  the  blade  row  exit,  because  of  the  low 
adiabatic  wall  temperatures  associated  with  high  Mach 
numbers.  At  high  power  levels,  and  thus  at  high  sur¬ 
face  temperatures,  the  adhesive  ability  of  the  double- 
faced  tape  is  reduced  and  the  heaters  could  unbond 
from  the  surface.  The  choice  of  crystal  temperatures, 
37.9®C  (100.2"F),  or  64.5"C  (148.1®F),  was  made  with 
consideration  to  these  high  and  low  power  limits. 

Definition  of  the  Heat  Transfer  Coefficient 
While  the  liquid  crystal  technique  offers  the  advantage 
of  good  spatial  resolution,  this  technique  does  not  fa¬ 
cilitate  heat  transfer  measurements  at  a  fixed  point  for 
varying  power  levels.  Measurements  of  this  type  can  be 
obtained  using  surface  mounted  thermocouples.  In  the 
present  study,  two  thin  film  thermocouples  were  bonded 
to  the  Inconel  heaters,  one  upstream  of  the  blade  row 
and  one  downstream  (see  Fig.  2).  Measurements  from 
these  thermocouples  were  used  to  determine  the  most 
appropriate  definition  of  heat  transfer  coefficient.  The 
thin  layers  of  paint  and  liquid  crystals  influenced  the 
thermocouple  measurements,  but  this  data  was  used 
only  to  compare  Stanton  number  definitions.  Adding 
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corrections  for  the  paint  conductance  simply  shifts  the 
level  of  but  the  trends  that  will  now  be  described 
remain  the  same.  Even  for  relatively  small  tempera¬ 
ture  differences  such  as  in  these  tests,  the  heat  trans¬ 
fer  coefficient  should  be  independent  of  the  heat  flux 
level.  Because  of  the  low  temperature  differences  and 
the  high  Mach  numbers,  the  local  Eckert  numbers  were 
very  high.  Consequently,  the  choice  of  temperature  dif¬ 
ference  used  as  the  driving  potential  has  a  significant 
effect  on  the  heat  transfer  coefficient.  The  most  appro¬ 
priate  definition  of  Stanton  number  is  the  one  which 
results  in  values  that  do  not  depend  on  the  level  of  sur¬ 
face  heat  flux.  In  general,  the  definitions  of  the  heat 
transfer  coefficient,  h,  and  of  the  Stanton  number,  St, 
are; 

a  _  q/A  _  h 
”  pUin^>^  “  pUinCp 

Each  of  the  terms  in  this  equation  can  be  defined 
or  evaluated  at  different  locations  or  at  different  tem¬ 
peratures.  In  the  following  discussion,  the  downstream 
thermocouple  readings  were  examined  to  determine  an 
appropriate  set  of  terms  for  use  in  calculating  the  Stan¬ 
ton  number.  These,  results  are  for  the  case  of  an  in¬ 
let  Rec^  =  0.973  x  10®,  Mex  =  1.328,  and  with  the 
turbulence  grid  installed.  This  case  provides  the  most 
contrast  between  definitions,  but  the  trends  are  consis¬ 
tent  for  both  thermocouples  and  for  all  flow  cases.  Five 
different  definitions  of  Stanton  number  are  described  in 
Table  2,  and  the  results  of  those  definitions  are  shown  in 
Fig.  4.  For  all  definitions  the  reference  velocity  was  cho¬ 
sen  as  Uin  and  Cp  was  assumed  to  be  constant.  Also  for 
all  of  the  definitions,  q  represents  the  net  surface  heat 
transfer  rate.  The  density  used  in  the  Stanton  num¬ 
ber  calculations,  unlike  the  value  used  for  the  Reynolds 
number,  is  not  the  actual  physical  density.  For  defini¬ 
tions  A  and  5,  p  is  the  inlet  total  density.  For  other 
definitions  it  is  calculated  using  inlet  total  pressure  and 
a  temperature  as  shown  in  Table  2.  The  local  adiabatic 
wall  temperature,  Taw^  is  defined  as: 

Tqw  _  . _ I  ^ _ 

1^“’’  1  +  0.5(7-1)M2 

The  isentropic  Mach  number,  M,  was  determined  from 
endwall  static  pressure  measurements.  A  reference  tem¬ 
perature  for  properties  given  by  Eckert  (1955)  also  ap¬ 
pears  in  Table  2.  Figure  4  clearly  shows  that  either  def¬ 
inition  J?  or  J?  is  the  most  appropriate  for  the  current 
study.  Both  of  the  definitions  unfortunately  require  an 
iterative  calculation  because  of  the  interdependency  be¬ 
tween  Taw  a-iid  Tref  •  Definition  E  accounts  for  the  bulk 
temperature  rise  caused  by  the  endwall  heaters  and  in 
general  would  therefore  be  more  appropriate  than  D. 


For  the  relatively  high  Reynolds  numbers  and  small 
temperature  differences  of  the  current  study,  bulk  heat¬ 
ing  has  little  effect  on  St,  but  at  lower  Reynolds  num¬ 
bers  or  larger  temperature  differences  this  may  not  be 
the  case.  Accounting  for  the  effects  of  bulk  heating  re¬ 
quires  an  integration  of  the  form  f  qdA  which  is  quite 
simple  for  a  uniform  q  boundary  condition,  but  which 
is  more  complex  if  the  flow  is  modeled  with  a  uni¬ 
form  wall  temperature  boundary  condition.  Since  many 
CFD  codes  implement  the  simpler  uniform  wall  temper¬ 
ature  boundary  condition,  all  data  presented  here  will 
be  based  on  definition  J9. 

Figure  4  also  shows  an  overall  variation  in  St  of 
about  5%  for  either  definition  D  or  E.  As  will  be  de¬ 
scribed,  the  overall  uncertainty  of  the  liquid  crystal 
measurements  was  estimated  to  be  6%  to  11%.  The  pri¬ 
mary  reason  for  the  high  sensitivity  of  the  Stanton  num¬ 
ber  to  the  choice  of  temperatures  is  that  local  Eckert 
numbers  are  high.  At  the  blade  row  exit,  and  at  a  tem¬ 
perature  corresponding  to  the  lower  temperature  liquid 
crystal,  the  loced  Eckert  number  was  approximately  8. 
At  these  relatively  high  Eckert  numbers,  heat  transfer 
rates  are  a  function  of  both  Reynolds  number  and  Mach 
number.  The  careful  definition  of  Stanton  number  de¬ 
scribed  above  properly  accounts  for  the  Mach  number 
dependency  through  the  use  of  the  adiabatic  wall  tem¬ 
perature  and  Eckert’s  reference  temperature.  The  refer¬ 
ence  temperature  accounts  for  thermophysical  property 
variations.  It  is  interesting  to  note  that  the  thermocou¬ 
ple  temperature  matched  T/c,!  at  qhtr  =  560  Watts  and 
matched  Tjc,2  at  qhtr  =  1300  Watts. 

The  net  surface  heat  transfer  rate,  q,  used  to  calcu¬ 
late  the  Stanton  number,  is  the  heater  power  corrected 
for  back-face  conduction  losses  and  for  radiative  heat 
transfer  to  other  surfaces: 

q  —  qhtr  f^foamA(Ts  Tfoam^/tfoam  €cA(T^ 

where  qhtr  =  VI  and  A  is  the  total  heater  area 
(824.2  cm^).  The  second  term  in  this  equation  rep¬ 
resents  the  back-face  heat  conduction  loss,  and  was  less 
than  0.6%  of  qhtr-  The  third  term  represents  the  radia¬ 
tive  losses,  and  was  less  than  4.5%  of  qhtr-  For  the  re¬ 
sults  shown  in  Fig.  4,  T,  was  measured  with  the  surface 
thermocouple.  When  endwall  Stanton  numbers  were 
being  determined,  Tg  was  the  temperature  of  the  liquid 
crystal.  Tjoam  was  measured  on  the  back  face  of  the  in¬ 
sulating  foam,  kfoam  and  tfoam  were  given  above,  and 
€  =  0.98. 
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Table  2  Definitions  used  for  Stanton 
number  terms 


qt  =  _ 

pUinCp{T.-Tg) 


Def’n 

P 

Tg  definition 

A 

P- 

Tf 

B 

p'i 

T  “  T/ 

T  = 

1  —  r 

l'’"''  1  + 0.5(7 -1)M2J 
Pri/3,  Pr  =  Pt{T,^) 

C 

PL 

Taw  —  T^n 

Trej 

r  = 

1  -  r  1 

RTrcf 

L'  '  1  +  0.5(7 -1)M2J 

=  r+o.5(r,-T)+ 
0.22(ra„  -  T) 
Pr^!^,  Pr  =  Pr{Tex) 

D 

P- 

Ta^  =  T!„ 

Trej 

r  = 

1  —  r 

L'  '  1  +  0.5(7  -  1)M2J 
=  r+o.5(r, -r)+ 
0.22(r<.„,  -  T) 
Pri/3,  pr  =  Pr(Tref) 

E 

P- 

. — 

ref 

Taw  —  C^in  “b  ^Thulk)^ 

1  + 0.5(7 -1)M21 
rre/  =  r+o.5(T,-r)+ 
0.22(r„„  -  T) 

r  =  Pr^!^,  Pr  =  Pr{Tre}) 

A  rr>  ^htrip^/ Cx  "b  0.3)T/i<r/A 

Up 

q  (Watts) 


note:  %  variations  defined  as  ^ —  x  100% 

Fig.  4  A  Comparison  of  Stanton 
number  definitions 


An  uncertainty  analysis  was  performed  following 
the  method  of  Kline  and  McClintock  (1953).  The  anal- 
ysis  showed  that  the  AT  term  in  the  calculation  of 
Stanton  number  was  the  major  source  of  uncertainty. 
Based  on  the  following  uncertainties:  T/c.i  =  ±0.5®C 
(±0.9"F),  Tic, 2  =  ±0.8"C  (±1.4"F),  and  Taw  =  ±0.6"C 
(±1.0®F),  the  overall  uncertainty  in  St  was  determined 
to  be  less  than  11%.  The  liquid  crystal  uncertainties 
were  estimated  by  applying  power  to  the  heaters  and 
monitoring  the  surface  mounted  thermocouples.  The 
crystals  were  viewed  using  the  same  photographic  tech¬ 
niques  that  were  used  during  actual  data  acquisition. 
To  minimize  conduction  errors  through  the  thin  paint 
layer  between  the  thermocouple  and  the  liquid  crys¬ 
tals,  air  flow  was  turned  off  and  the  pressure  in  the  test 
section  was  reduced  to  near  vacuum,  thus  minimizing 
natural  convection.  During  the  actual  data  acquisition, 
the  high  temperature  liquid  crystal  data  were  used  if 
the  temperature  difference,  AT,  using  the  low  temper¬ 
ature  liquid  crystal  was  too  low  to  maintain  the  desired 
maximum  uncertainty  in  Stanton  number. 

MEASUREMENT  RESULTS 
Inlet  Flow  Measurements  Well-documented  inlet 
boundary  conditions  are  needed  for  CFD  calculations. 
Figure  5  shows  the  results  of  aerodynamic  inlet  mea¬ 
surements  on  a  plane  one  axial  chord  upstream  of  the 
blade  leading  edge  plane.  Figures  5a  and  5b  show  local 
Mach  number  contours  without  a  turbulence  grid  and 
with  the  turbulence  grid  installed.  Both  figures  present 
data  obtained  at  Rec^  =  1.0  x  10®.  Based  on  flat  plate 
boundary  layer  correlations,  the  Rec^  =  0.5  x  10®  cases 
can  approximately  be  expected  to  have  a  10%  thicker 
inlet  boundary  layer.  The  exit  Mach  number  for  the 
end  wall  heat  transfer  cases  was  either  1.0  or  1.3  leaving 
the  flow  choked  at  the  throat  and  thus  leaving  the  inlet 
flow  independent  of  M^x  • 

Figure  5  shows  that  without  the  turbulence  grid 
there  is  a  thick  and  pitchwise  periodic  inlet  boundary 
layer.  Because  of  the  large  leading  edge  diameter,  the 
effects  of  the  blades  are  evident  even  one  axial  chord 
upstream  of  the  leading  edge.  Comparing  Figs.  5a  and 
5b  shows  that  the  presence  of  the  grid  adversely  af¬ 
fects  the  flow  uniformity.  The  primary  measurement 
passage,  Passage  5,  has  a  midrange  level  of  turbulence 
intensity,  but  a  higher  than  average  velocity.  In  terms 
of  Reynolds  number,  Rec,.  in  Passage  6  was  reduced 
about  8%  while  Rec^  in  Passage  4  remained  about 
the  same.  Without  a  grid,  the  inlet  Tu  was  uniform 
at  about  0.25%  for  Rec^  =  1.0  x  10®  and  0.50%  at 
Rec^  =  0.5  X  10®.  Figure  5c  shows  the  turbulence  in¬ 
tensity  distribution  with  the  grid  installed.  The  turbu- 
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lence  intensity,  Ttx,  in  Fig.  5c  has  been  normalized  by 
the  spanwise  and  pitchwise  averaged  inlet  velocity.  The 
roughly  linear  dependence  of  Tu  on  pitchwise  distance 
resulted  from  the  grid  being  installed  perpendicular  to 
the  mean  flow,  thus  making  it  not  parallel  to  the  mea¬ 
surement  plane  and  the  blade  leading  edge  plane.  The 
measurement  and  blade  leading  edge  planes  were  par¬ 
allel.  Gregory-Smith  and  Cleak  (1992)  showed  that  a 
42.75®  skewed  grid  produced  a  mean  flow  deflection  of 
more  than  1.0®.  Skewing  the  grid  in  the  present  study 
may  have  resulted  in  a  more  uniform  Tu  distribution, 
but  it  probably  would  have  had  an  even  more  detri¬ 
mental  effect  on  the  mean  velocity.  The  spanwise  and 
pitchwise  averaged  inlet  velocity  was  used  to  calculate 
Rec^  and  St, 


05  q 


0  -1  y/span  -2  -3 

a,  Mach  number  distribution  without  turbulence  grid 


b,  Mach  number  distribution  with  turbulence  grid 


y/span 

c,  turbulence  intensity  with  turbulence  grid 


Fig.  5  Aerodynamic  inlet  conditions 

Turbulence  length  scales  can  have  a  significant  in¬ 
fluence  on  heat  transfer  rates.  A  longitudinal,  integral 
turbulence  length  scale  for  the  present  case  was  esti¬ 
mated  by  extrapolating  data  given  by  Van  Fossen  et  al. 
(1994)  to  higher  air  speeds.  Their  data  was  taken  in 
a  low  speed  wind  tunnel  (45  m/s)  of  the  same  width 
(15.24  cm)  and  with  the  same  grid  configuration  (bar 
size  and  spacing)  as  the  present  test.  As  shown  by 
Roach  (1987)  and  verified  by  Van  Fossen  et  al.  (1994), 


the  length  scale  generated  by  grids  is  independent  of 
Reynolds  number,  therefore,  extrapolation  to  higher 
speeds  is  reasonable.  The  length  scale,  LxfCx,  thus 
estimated  was  0.76  at  the  leading  edge  plane  of  Pas¬ 
sage  5. 

Endwall  Static  Pressure  Measurements  As  previ¬ 
ously  discussed,  static  pressure  endwall  distributions 
were  needed  along  with  the  surface  temperatures  to 
evaluate  Stanton  numbers.  Figure  6  shows  the  static-to- 
inlet  total  pressure  ratio  measurements  for  an  exit  Mach 
number  of  1.3.  Also  shown  are  predictions  made  us¬ 
ing  the  three-dimensional  Navier-Stokes  code,  RVC3D, 
(Chima  and  Yokota,  1990).  The  analysis  was  done  for 
flow  around  a  single  blade.  The  graph  of  the  endwall 
pressure  distribution  was  then  offset  by  one  pitch  and 
repeated  for  comparison  with  data.  Excellent  agree¬ 
ment  is  seen  between  the  data  and  the  calculations 
with  the  only  substantial  difference  being  with  the  0.3 
contours  downstream  of  the  trailing  edge.  The  mini¬ 
mum  P/Pin  values  for  the  data  and  for  the  calculations 
were  0.3055  and  0.2796  respectively.  These  values  cor¬ 
respond  to  only  a  0.6®C  (1®F)  difference  in  Taw^  The 
dots  in  the  experimental  data  figure  show  the  location 
of  the  pressure  taps.  Even  with  87  taps  per  passage, 
the  experimental  data  may  not  have  enough  resolution 
to  adequately  define  these  low  pressure  islands.  For 
this  reason,  the  isentropic  Mach  number  used  to  calcu¬ 
late  the  Stanton  number  was  obtained  from  the  CFD 
calculations  rather  than  directly  from  the  data.  The 
calculated  Mach  number  distribution  was  not  affected 
by  a  Reynolds  number  vaxiation,  so  only  two  exit  Mach 
number  calculations  were  needed  to  process  the  data. 
Other  endwall  static  pressure  data  is  available  from  Giel 
et  al.  (1996). 


a  experimental  data  with  tap  locations 


b.  CFD  calculations 

Fig.  6  Endwall  static  pressure  distributions 
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Figure  6  also  shows  the  strong  three-dimension¬ 
ality  of  the  flow.  This  is  particularly  evidenced  by  the 
distortion  of  the  P/P/yj  =  0.8  contour  line.  The  distor¬ 
tion  is  caused  by  the  horseshoe  vortex.  Pressure  mea¬ 
surements  on  the  blade  surface,  Giel  et  al.  (1996),  also 
showed  strong  three-dimensionality  in  the  flow  field. 
Near  midchord  there  was  nearly  a  thirty  percent  de¬ 
crease  in  the  rotor  pressure  differential  at  6%  of  span 
compared  to  the  midspan  value. 

Heat  Transfer  Measurements  Heat  transfer  mea¬ 
surements  were  made  at  two  Reynolds  numbers  and 
at  each  of  two  exit  Mach  numbers  with  a  low  turbu¬ 
lence  inlet  condition.  Each  of  these  tests  were  repeated 
using  the  grid  to  generate  a  high  level  of  turbulence. 
Results  are  presented  in  this  paper  for  the  full  facto¬ 
rial  experiment  of  three  variables  at  two  levels  each, 
i.e.,  eight  cases.  Table  3  lists  test  conditions  for  the 
eight  cases.  Note  that  one  of  the  low  Reynolds  number 
cases,  Case  5,  showed  a  greater  variability  in  exit  con¬ 
ditions  during  the  run.  However,  even  this  variability 
was  small. 

Table  3  Description  of  endwall  heat 
transfer  cases 


Case 

Mex 

Tu  grid 

1 

0.937  ±  0.007  X  lO” 

1.321  ±  0.002 

no 

2 

0.938  ±  0.009  X  10® 

0.962  ±0.001 

no 

3 

0.973  ±  0.016  X  10® 

1.328  ±0.001 

yes 

4 

0.952  ±  0.017  X  10® 

0.992  ±  0.001 

yes 

5 

0.473  ±  0.004  X  10® 

1.272  ±0.036 

no 

6 

0.475  ±  0.003  X  10® 

0.961  ±0.004 

no 

7 

0.493  ±  0.003  X  10® 

1.313  ±0.006 

yes 

8 

0.493  ±  0.004  X  10® 

0.982  ±0.002 

yes 

{all  repeatabilities  based  on  95%  confidence  limits) 


The  endwall  Stanton  number  distributions  for  the 
eight  cases  are  presented  in  Figs.  7  through  14.  Each 
of  these  figures  will  be  first  discussed  separately,  then 
compared  to  each  other  in  order  to  summarize  the  ef¬ 
fects  of  Reynolds  number,  exit  Mach  number,  and  the 
inlet  turbulence  grid.  The  figures  are  arranged  such 
that  the  effects  of  Reynolds  number  can  be  observed  by 
comparing  any  two  figures  that  are  horizontally  oppo¬ 
site  each  other,  e.g.,  Fig.  7  and  Fig.  11.  The  effects  of 
Mach  number  can  be  observed  by  comparing  sequen¬ 
tially  numbered  figures,  e.g.,  Fig.  7  and  Fig.  8.  Finally, 
the  effects  of  the  turbulence  grid  can  be  observed  by 
comparing  every  other  figure,  e.g.,  Fig.  7  and  Fig.  9. 
Because  all  of  the  data  was  interpolated  onto  the  same 
rectangular  grid,  it  was  possible  to  subtract  the  local 
values  from  any  two  cases  in  order  to  determine  per¬ 
centage  changes  in  St  brought  about  by  changes  in  one 


or  more  flow  parameters.  While  the  authors  used  this 
procedure  to  help  interpret  the  results,  the  large  num¬ 
ber  of  figures  generated  from  this  procedure  prevents 
their  inclusion  here.  Simply  subtracting  St  values  at 
the  same  {x,y)  location  from  any  two  figures  leads  to 
the  same  conclusions. 

Figure  7  shows  the  endwall  Stanton  number  dis¬ 
tribution  for  a  case  with  an  inlet  Rec^.  =  1.0  x  10^,  a 
PexIPin  corresponding  to  an  isentropic  exit  Mach  num¬ 
ber  of  1.3,  and  without  the  turbulence  grid  -  Case  1. 
The  heated  section  begins  0.3  upstream  of  the  blade 
leading  edge  plane.  High  Stanton  number  levels  are  ex¬ 
pected  in  this  region  because  of  the  unheated  starting 
length  heat  transfer  effects.  Even  though  the  Stanton 
number  is  in  theory  infinite  at  the  heater  leading  edge,  a 
flat  plate  correlation  (Kays  and  Crawford,  1980)  shows 
that  it  is  reduced  to  2.4  x  10“^  at  a  distance  O.lCx 
from  the  start  of  heating.  The  low  value  at  this  loca¬ 
tion  results  from  the  thick  inlet  boundary  layer  which 
corresponds  to  a  Reynolds  number,  of  14.6  x  10®. 

Further  downstream,  the  presence  of  the  blades 
quickly  distorts  the  Stanton  number  distribution  from 
the  flat  plate  correlation  results.  Near  the  blade  lead¬ 
ing  edge,  the  Stanton  numbers  are  significantly  influ¬ 
enced  by  the  horseshoe  vortex  structure.  As  the  flow 
approaches  the  leading  edge,  there  is  a  general  increase 
in  heat  transfer.  High  Stanton  number  levels  are  seen 
very  close  to  the  blade,  near  its  stagnation  point,  in 
agreement  with  the  effects  observed  by  Goldstein  and 
Spores  (1988)  of  a  leading  edge  corner  vortex.  As  the 
flow  accelerates  around  the  suction  surface,  there  is  a 
region  of  very  high  heat  transfer,  St  x  1000  >  4.5.  In 
the  pitchwise  direction,  this  is  a  region  of  high  gradi¬ 
ents,  since  near  midpassage,  there  is  a  relatively  low 
heat  transfer  region. 

A  comparison  of  the  data  in  the  primary  measure¬ 
ment  passage  with  the  data  in  the  adjacent  passage  in¬ 
dicates  that  the  data  shows  good  periodicity.  Minimum 
heat  transfer  regions  are  seen  at  midchord  near  each 
of  the  pressure  surfaces.  This  is  as  expected,  since  as 
Fig.  6  shows,  this  is  a  low  velocity  region.  From  here, 
Stanton  number  levels  gradually  increase  going  down¬ 
stream,  Near  the  suction  surface,  levels  decrease  in  the 
streamwise  direction  from  the  peak  levels  downstream 
of  the  stagnation  region.  The  effects  of  the  pressure- 
side  leg  of  the  horseshoe  vortex  are  seen  clearly  further 
along  the  suction  surface.  When  the  pressure-side  leg  of 
the  horseshoe  vortex  lifts  off  of  the  endwall,  it  no  longer 
scours  the  endwall.  A  small  peak  in  St  of  4.5  x 
is  observed  near  the  presumed  lift-off  point.  From  this 
point,  the  passage  vortex  causes  an  accumulation  of 
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secondary-flow  fluid  near  the  suction  surface.  This  fluid 
has  passed  over  the  heated  endwall,  thus  resulting  in 
lower  heat  transfer  rates  near  the  suction  surface.  In 
fact,  the  lowest  Stanton  number  levels  are  seen  on  the 
endwall  downstream  of  the  throat  along  the  uncovered 
portion  of  the  suction  surface.  Note  that  the  minimum 
Stanton  number  in  this  region  is  somewhat  lower  in  the 
primary  measurement  passage  (Passage  5)  than  in  the 
adjacent  passage  (Passage  6).  This  difference  is  caused 
by  the  fact  that  heaters  cover  the  entire  secondary  flow 
path  in  Passage  5,  but  only  a  portion  of  Passage  6  is 
covered  by  heaters.  This  allows  the  secondary  flow  to 
become  warmer  in  Passage  5  before  it  reaches  the  suc¬ 
tion  surface,  thus  decreasing  the  rate  of  heat  transfer. 

Highly  localized  peak  endwall  heat  transfer  rates 
occur  in  the  wake  region,  0.20  to  0.25  Ca;  downstream 
of  the  blade  trailing  edge.  This  distance  corresponds 
to  approximately  5  to  6  trailing  edge  diameters.  Gold¬ 
stein  and  Spores  (1988)  saw  a  similar  highly  localized 
peak,  but  for  their  tests  it  was  much  closer  to  the  blade, 
only  0.06  to  0.07  Cx  downstream  of  the  blade  trailing 
edge.  Further  downstream  of  the  blade  row  several 
distinct  regions  of  low  heat  transfer  are  seen.  These 
regions  are  probably  defined  by  the  shock  structure 
evident  in  Fig.  6.  The  Stanton  number  generally  de¬ 
creases  in  the  streamwise  direction  downstream  of  the 
blade  row  as  the  thermal  boundary  layer  grows  on  the 
endwall.  In  general,  Fig.  7  shows  good  endwall  heat 
transfer  periodicity  and  strong  evidence  of  the  highly 
three-dimensional  flow  field. 

The  endwall  Stanton  number  distribution  for 
Case  2  is  shown  in  Fig.  8.  The  conditions  for  this  case 
were  the  same  as  for  Case  1  except  that  the  isentropic 
exit  Mach  number  was  reduced  to  slightly  less  than  1.0 
so  that  no  shocks  would  be  present  in  the  downstream 
region.  The  mass  flows  for  Case  1  and  Case  2  were 
nearly  identical.  The  endwall  heat  transfer  levels  up¬ 
stream  of  the  throat  were  expected  to  be  the  same  as 
for  Case  1.  Comparing  Figs.  7  and  8  shows  that  this 
is  generally  true  within  the  experimental  uncertainty. 
The  figure  does  not  quite  show  the  same  peak  St  of 
4.5  X  10“^  near  the  blade  stagnation  point  as  was  seen 
in  Case  1.  This  was  probably  due  to  insufficient  resolu¬ 
tion  of  qhir  settings  when  acquiring  the  data  for  Case  2. 
Along  the  pressure  surface,  St  increases  gradually  in 
streamwise  direction  with  contour  lines  approximately 
being  normal  to  the  pressure  surface.  This  was  also  true 
for  Case  1. 

Because  there  are  no  shocks  present  to  interact  with 
the  endwall  boundary  layers,  Stanton  number  contours 
from  midchord  to  the  exit  generally  follow  the  paths 


of  expected  inviscid  streamlines  and  do  not  appear  in 
distinct  regions  as  they  did  in  Case  1.  The  low  heat 
transfer  region  along  the  uncovered  portion  of  the  suc¬ 
tion  surface  is  still  present  but  with  higher  St  levels 
and  with  a  smaller  area.  The  level  of  the  peak  heat 
transfer  {St  >  5.5  x  10“^)  downstream  of  the  trailing 
edge  remains  the  same,  but  its  location  is  slightly  closer 
to  the  trailing  edge,  approximately  3  to  5  trailing  edge 
diameters  downstream.  Downstream  of  the  blade  row 
the  Stanton  number  levels  are  40%  to  60%  higher  than 
Case  1  and  generally  follow  the  blade  wakes  without 
being  perturbed  by  shocks. 

Figure  9  shows  the  endwall  Stanton  number  dis¬ 
tribution  for  the  same  Reynolds  number  and  exit  Mach 
number  as  Fig.  7  but  obtained  with  the  inlet  turbulence 
grid  installed.  In  general,  the  figure  shows,  somewhat 
surprisingly,  that  the  primary  effect  of  the  turbulence 
grid  is  to  lower  peak  heat  transfer  levels.  The  most 
significant  differences  are  in  the  regions  in  front  of  the 
blade  leading  edges,  in  the  path  of  the  horseshoe  vor¬ 
tex.  The  low  Stanton  number  regions  of  2.5  x  10”^  are 
still  seen  upstream  of  this  region.  The  steep  St  gradi¬ 
ents  near  midpitch  seen  in  Fig.  7,  however,  are  missing. 
Recall  from  Fig.  5  that  the  turbulence  grid  significantly 
reduced  the  inlet  boundary  layer  thickness.  Blair  (1983) 
showed  that,  when  corrected  for  the  length  scale  a 
freestream  turbulence  intensity  of  7%  can  be  expected 
to  result  in  less  than  a  10%  increase  in  St.  The  reduced 
secondary  flows  resulting  from  the  thinner  inlet  bound¬ 
ary  layer  may  overshadow  this  level  of  increased  St.  The 
peak  St  area  seen  in  Fig.  7  near  the  presumed  vortex 
lift-off  point  has  moved  away  from  the  suction  surface, 
followed  by  10%  to  20%  higher  Stanton  numbers  along 
the  remainder  of  the  suction  surface.  This  observation 
supports  the  idea  that  the  secondary  flow  structure  is 
different  for  the  two  cases.  Figure  9  shows  that  the  re¬ 
gion  downstream  of  the  blade  trailing  edge  is  similar 
for  the  two  cases.  This  is  probably  due  to  the  fact  that 
even  for  Case  1,  the  horseshoe  vortex  structure  scours 
the  endwall,  increasing  the  freestream  turbulence  levels 
and  eliminating  the  thick  inlet  boundary  layers.  Down¬ 
stream  of  the  horseshoe  vortex  therefore,  the  two  cases 
should  have  similar  endwall  heat  transfer  distributions. 
The  far  downstream  region,  in  fact,  supports  this  sup¬ 
position  with  levels  and  patterns  very  similar  to  those 
of  Case  1. 

Figure  10  shows  the  final  high  Reynolds  number 
case,  Case  4.  The  trends  and  differences  for  this  case 
are  generally  the  same  as  for  the  cases  discussed  above. 
In  comparing  Fig.  10  with  Fig.  8,  it  is  evident  that  the 
increased  St  levels  near  the  suction  surface  are  more 
pronounced  for  this  case.  Downstream  of  the  blade  row, 
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however,  the  Stanton  number  levels  are  again  very  sim« 
ilar. 

Figure  11  shows  the  endwall  Stanton  number  dis¬ 
tribution  for  Case  5.  Contrasting  this  figure  with  Fig.  7 
shows  the  effects  of  reducing  the  Reynolds  number  for 
a  supersonic  exit  case  with  no  inlet  turbulence  grid. 
These  effects  are  significant  in  the  upstream  and  lead¬ 
ing  edge  region.  Increases  in  Si  are  evident  in  excess  of 
the  15  to  20%  increase  that  would  be  expected  assum¬ 
ing  that  St  ^  based  on  a  turbulent  flow  anal¬ 

ogy.  The  peak  Stanton  number  in  the  stagnation  region 
increases  from  about  4.5  x  10“^  for  Case  1  to  about 
5.5  X  10“*^  for  Case  4.  Aside  from  the  peak  value,  the 
overall  shape  of  the  contours  in  the  leading  edge  region 
remains  about  the  same,  resulting  in  higher  gradients 
of  St  near  midpitch.  The  blade  passage,  particularly 
near  the  pressure  surface,  shows  the  most  significant 
Reynolds  number  effects.  The  lower  Reynolds  number 
moves  the  minimum  St  region  significantly  closer  to 
the  suction  surface  and  slightly  downstream.  Stanton 
number  levels  at  quarter-chord  near  the  pressure  sur¬ 
face  increase  by  more  than  50%.  Levels  along  the  suc¬ 
tion  surface  also  increase  by  10  to  30%.  However,  levels 
near  midpassage,  from  0.5  to  0.75  Cx,  decrease  by  30% 
or  more.  This  all  has  the  effect  of  producing  very  large 
St  gradients  near  midchord.  These  large  effects  are  al¬ 
most  certainly  due  to  the  thicker  inlet  boundary  layer 
and  the  changes  in  the  horseshoe  vortex  that  it  brings 
about.  Downstream  of  the  throat,  the  uncovered  suc¬ 
tion  surface  region  and,  in  fact,  the  entire  downstream 
region  appears  to  show  levels  that  are  generally  higher 
than  for  the  Rec^  =  1.0  x  10®  case.  The  peak  levels 
downstream  of  the  trailing  edge  are  again  nearly  the 
same,  but  the  area  of  higher  heat  transfer  surrounding 
this  peak  is  significantly  larger. 

Figure  12  shows  the  endwall  Stanton  number  dis¬ 
tribution  for  Case  6.  The  area  upstream  of  the  throat 
agrees  reasonably  well  with  the  data  of  Case  5.  The  area 
downstream  of  the  throat  generally  shows  increased  St 
levels,  the  same  M^x  effect  that  was  seen  for  other  cases. 
In  comparing  Fig.  12  with  Fig.  8,  it  is  seen  that  the 
Reynolds  number  effects  axe  also  consistent  with  those 
discussed  above. 

The  endwall  Stanton  number  distribution  for 
Case  7  is  shown  in  Figure  13.  This  is  a  0,5  x  10® 
Reynolds  number  case  with  the  inlet  turbulence  grid 
installed.  The  lower  Reynolds  number  would  normally 
produce  a  thicker  inlet  boundary  layer,  but  the  tur¬ 
bulence  grid  reduces  it  significantly.  The  net  result  of 
these  effects  is  that  the  low  pressure  surface  St  region 
characteristic  of  the  higher  Reynolds  number  cases  has 


still  moved  downstream  and  closer  to  the  suction  sur¬ 
face,  but  not  as  much  as  for  the  no-grid  cases  (5  and 
6).  Therefore,  the  gradients  of  St  near  the  suction  sur¬ 
face  at  midchord  are  not  as  large  as  for  the  no-grid 
cases.  Decreasing  the  Reynolds  number  with  the  tur¬ 
bulence  grid  installed  produces  trends  similar  to  the 
no-grid  cases,  but  the  magnitude  of  the  differences  is 
reduced. 

The  final  endwall  Stanton  number  distribution, 
Case  8,  is  seen  in  Fig.  14.  Comparing  Fig.  14  to  Fig.  10 
shows  the  same  Reynolds  number  trends  that  have  been 
shown  throughout.  Comparing  it  to  Fig.  12  shows 
consistent  turbulence-grid  effects.  When  compared  to 
Fig.  13,  the  consistent  downstream  Mach  number  trend 
of  40%  to  60%  higher  St  levels  associated  with  the  lower 
Mex  levels  is  again  seen.  Also,  the  agreement  with 
Case  7  in  the  region  upstream  of  the  throat  is  excel¬ 
lent. 

SUMMARY  AND  CONCLUSIONS 

In  general,  lowering  iiec*  by  a  factor  of  two  has  sig¬ 
nificant  effects  on  the  endwatll  heat  transfer  in  excess  of 
the  15  -  20%  increases  expected  from  flat  plate  correla¬ 
tions.  Stanton  number  levels  near  the  pressure  surface 
increase  by  more  than  50%,  while  levels  near  the  suction 
surface  increase  by  nearly  30%.  Lower  Reynolds  num¬ 
bers  also  move  the  minimum  passage  Stanton  number 
slightly  downstream  and  towards  the  suction  surface, 
decreasing  levels  at  midpassage  by  30%  or  more.  These 
effects  were  shown  to  result  in  large  St  gradients  near 
the  suction  surface  at  about  midchord.  The  Reynolds 
number  effects  were  most  evident  in  Cases  5  and  6  with 
no  grid,  but  were  also  seen  for  Cases  7  and  8  with  the 
turbulence  grid. 

The  shock  structure  present  in  the  supersonic  M^x 
cases  were  shown  to  break  up  the  smooth  St  distribu¬ 
tions  present  in  the  sonic  Mex  cases.  In  general,  Stanton 
number  levels  are  40%  to  60%  higher  for  the  lower  exit 
Mach  number  cases.  In  all  cases,  the  downstream  Stan¬ 
ton  number  levels  decrease  in  the  streamwise  direction 
as  the  thermal  boundary  layer  grows  on  the  endwall. 

In  general,  the  data  show  the  somewhat  surprising 
result  that  the  primary  effect  of  the  turbulence  grid  is 
to  lower  peak  heat  transfer  levels.  This  result  is  thought 
to  be  due  to  the  secondary  flow  change  brought  about 
by  the  reduced  inlet  boundary  layer  thickness.  Down¬ 
stream  of  the  blade  row  the  turbulence  grid  had  essen¬ 
tially  no  effect  on  the  endwall  heat  transfer. 

Overall,  the  steady-state  liquid  crystal  technique 
proved  to  be  quite  valuable  for  obtaining  CFD  code 
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validation  endwall  heat  trzoisfer  data.  The  analysis  of 
thermocouple  data  showed  the  importance  of  a  careful 
Stanton  number  definition.  The  good  spatial  resolu¬ 
tion  offered  by  the  large  scale  of  the  cascade  and  by 
the  liquid  crystal  technique  allowed  the  effects  of  the 
large  secondary  flows  in  the  blade  passages  to  be  clearly 
quantified.  The  data  show  the  significant  differences  in 
the  heat  transfer  distributions  that  resulted  from  varia¬ 
tions  in  Reynolds  number,  exit  Mach  number,  and  inlet 
turbulence  grid.  For  each  of  the  eight  cases,  however, 
a  peak  heat  transfer  region  occurred  in  a  small  down¬ 
stream  region,  between  3  to  6  diameters  downstream  of 
the  trailing  edge. 

Well-defined  inlet  flow  measurements  were  pre¬ 
sented,  important  for  specifying  boundary  conditions  to 
CFD  codes.  The  endwall  heat  transfer  data  presented 
here,  along  with  the  aerodynamic  data  presented  by 
Giel  et  al.  (1996)  comprise  a  complete  set  of  data  suit¬ 
able  for  CFD  code  and  model  validation.  Also,  as  re¬ 
quired  by  most  CFD  analyses,  the  endwall  heat  transfer 
measurements  show  good  periodicity.  Electronic  tabu¬ 
lations  of  all  eight  cases  presented  in  this  paper  are 
available  upon  request.  Comparing  CFD  calculations 
to  all  of  the  eight  test  conditions  available  would  give 
greater  confidence  in  the  ability  of  CFD  to  predict  end- 
wall  heat  transfer. 
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